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Effects of dynamical quarks on the microscopic spectrum of the Wilson Dirac operator are ana- 
lyzed by means of effective field theory. We consider the distributions of the real modes of the 
Wilson Dirac operator as well as the spectrum of the Hermitian Wilson Dirac operator, and work 
out the case of one flavor in all detail. In contrast to the quenched case, the theory has a mild 
sign problem that manifests itself by giving a spectral density that is not positive definite as the 
spectral gap closes. 
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1. Introduction 

The low-lying spectrum of the Dirac operator is a fascinating subject, which also contains 
important physics. Recently, three of the present authors presented an analysis of the quenched 
spectrum of the Wilson Dirac operator Dw Pp. Focus was on the low-lying real modes of Dw and 
the Hermitian counterpart D5 = 75 {Dw+m). Effects of the lattice spacing a were taken into account 
to leading order. A chiral Random Matrix Theory that encapsulates these leading order terms was 
then established. As the size N of the random matrices goes to infinity, a scaling regime is reached 
where this chiral Random Matrix Theory coincides with what one obtains from the Wilson chiral 
Lagrangian to leading order in a. This extends in a precise manner the universal Random Matrix 
Theory results of continuum fermions [Q, ||] to Wilson fermions in the microscopic scaling regime. 
In particular, the appropriate definition of an £-regime [Q] for the low-lying eigenvalues of the 
Wilson Dirac operator is identified [JTJ| _ The corresponding spectrum away from the microscopic 
limit was first analyzed at the mean field level by Sharpe in ref. [ph. The Letter §\§ was very much 
motivated by that work and a wish to understand in detail and at an analytical level some of the 
results of the lattice simulations in ref. [Q]. 

Here we report on a study of the effect of dynamical quarks on these results. Because the case 
of two light flavors is significantly more difficult in terms of computational complexity, we take 
here the first step of unquenching by considering Nf = 1 . This case is of interest in its own right 
because there are no Goldstone bosons and hence no chiral Lagrangian at our disposal. Never- 
theless, effective field theory can be used to describe in a precise way the leading-order effects of 
Wilson terms in lattice gauge theory also in this case. By projecting onto sectors of a fixed num- 
ber v of real modes (counted with the sign of their chiralities, see below), we can also establish a 
chiral Random Matrix Theory with exactly the same properties as the effective field theory in the 
scaling limit. As for continuum fermions, the effective field theory in each fixed sector looks just 
like the leading term in an e-regime counting of a chiral Lagrangian. Yet there are no Goldstone 
bosons and hence no way to systematically introduce a full-fledged space-time dependent chiral 
Lagrangian which could incorporate sub-leading effects of an associated £-expansion. 

2. The effective field theory 

Chiral symmetry for QCD with just one flavor is broken explicitly due to the U (1) anomaly, 
and there are no Goldstone bosons. As a consequence, we do not have the toolbox of chiral Pertur- 
bation Theory available. Leutwyler and Smilga [Q] faced a similar situation when dealing with the 
spectrum of the continuum Dirac operator, and we will here follow the same line of reasoning. In 
the continuum, the leading effect of a quark mass m is proportional to the four-volume V. Because 
the logarithmic derivative yields the chiral condensate E, it follows that the partition function must 
read Z ~ expfraiy]. This term corresponds to m\j/\j/ in the QCD Lagrangian. For Wilson fermions, 
the Symanzik effective action has additional operators ~ cr{\\f\\f) . Such terms give an additional 
contribution to the free energy of order a 2 so that now 

Z = exp [mZV -2W 8 Va 2 ] (2.1) 

where Wg is a so far unknown constant. We have chosen the parametrization so that this constant 
is naturally positive (the factor of 2 is for later convenience). As argued in ref. [|l]] a positive sign 
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of W& is follows from from the Hermiticity properties of Dy?. We define the appropriate £-regime 
here by requiring that 

m = rriLV and a 2 = a 2 W^V 

remain fixed as V — > °°. This is the regime where there is competition between m and a 2 effects 
and one can imagine that a phase transition may occur. This turns out to be a transition to the 
Aoki phase []7|] (see also [Q]). Other countings can also be considered [Q], but they are not of direct 
interest to us here. 

In the continuum, a chiral rotation a shifts the vacuum angle 8 — > 6 + a. Noting that the 
a 2 -term in the effective action comes from operators a 2 (\j/\j/) 2 , we define 

Z(6) = exp [mcos (d)LV -2W 8 Va 2 cos(2d)] , (2.2) 

and its Fourier components read: 

— e' ve exp [mcos(d)ZV - 2W s Va 2 cos(2d)] (2.3) 
-n 2k 

Inverting this, we recover the original partition function as a sum over each Z v after setting 6=0: 

Z(0 = O) = £ Z v . (2.4) 



Let us now look at each Z v separately. Interestingly, 



Z v = [ dU det(£/) v exp 

Jum 



-mIVTr[U + U- 1 } - W s Va 2 Tr[U 2 + U~ 2 } 



(2.5) 



This looks exactly like the zero-momentum piece of the leading terms of a U ( 1 ) chiral Lagrangian 
for Wilson fermions [|10[]. However, there are no Goldstone bosons, and the U(l) 'degree of free- 
dom' results from the angular integration variable of the Fourier transform. 

For general Nf there would also be double-trace terms like (Tr [U + U ~ 1 ] ) 2 and (Tr[U — U~ 1 ]) 2 , 
but in this U(l) case such terms just change the normalization of Wg, after use of elementary trigono- 
metric identities. 



3. Low-lying modes of the Wilson Dirac Operator 

To get spectral information for the Wilson Dirac operator we need either 

* Pairs of extra species with opposite statistics (the graded method Jl2|]) or 



* Replicas ||11 



Here we use the graded method. We thus add a bosonic quark and a corresponding additional 
fermionic quark, both with appropriate sources. When these sources are set equal to each other, the 
two additional determinants exactly cancel each other. In this limit, the partition function of this 
graded theory therefore equals the partition function of QCD with the original one flavor. 

The graded method can be used in the effective field theory as well. Additional Grassmann 
integrations truncate and trivially converge, but care must be taken to ensure convergence of the 
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bo sonic integrations. This problem has been solved in the context of continuum fermions in ref. 



Jl2[]. The graded partition function is 

dUSd & t(U) v e iL 2 ST <^ u - u ~ 



Z2u(^,jr) 



The source terms are 
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l ])+d 2 STt(U 2 +U- 2 ) 
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and when m = m' and z = z' a little miracle occurs: the graded partition function becomes equal 
to the original partition function of Nf = 1. This follows from general principles, but it arises in a 



highly non-trivial manner from the actual integrations of eq. (3.1). An explicit parametrization of 



the graded matrix U has been provided in ref. [ ]12| ] : 

? it+iu cos(6) ie it+i ^>sm(e 
U - I ie*-^sin(0) e^ itt cos(0 





exp 



/0 Ooi 
a 2 
\Pi ft 



where d,t,u G [ — tt, 7t] and (j> £ [0, 7t]. The bosonic degree of freedom 5 is integrated over the real 
line, and the a's and j8's are Grassmann variables. 



A careful reader will have noticed the unusual form of (3.1). Before extending the theory to 
the graded case, a rotation U — > iU has been performed. In the original U (l)-integral this simply 
shifts the angular variable by n/2, while still integrating it over the full circle. Doing such a rotation 
prior to extending the action to the graded case corresponds to a particular path of integration for the 
bosonic variable s. It is the integration path used in eq. (fO|) which corresponds to a non-Hermitian 
(but 75-Hermitian) Wilson Dirac operator D w . 



The integrals in eq. ( |3.1| ) are tedious but doable. We have performed the Grassmann integra- 
tions and one of the angular integrations explicitly. The resulting expressions will be published 
elsewhere. Here we choose to present our results in a graphical manner. 

While the Wilson Dirac Operator is not Hermitian, it is important that it nevertheless retains 
75-Hermiticity: D^y = Y5DwY5- Indeed, it is this property that ensures Hermiticity of D5. The 
spectrum of Dw thus lies in the complex plane, each non-real eigenvalue being matched by its 
complex conjugate partner. To compute the spectrum of the (non-Hermitian) Wilson Dirac operator 
by analytical means is slightly cumbersome because of this. However, Dw also has a certain number 
of eigenvalues sitting on the real line. The distribution of the chiralities of the corresponding states 
over the Dirac spectrum is much easier to compute. To this end, let us define a resolvent (and put 
z = t = 0) 

L v (rhf,m) = lim -^XnZZAmf^^m'). (3.2) 

A'-»h am 1 

The discontinuity across the real line gives us the distribution of the chiralities over the Dirac 
spectrum 

P V X CQ^ I 8{£-&)Zk = ^HL v (* f ,t)] (3.3) 



4 



Microscopic Wilson Dime spectrum for Nf = 1 



Poul H. Damgaard 




with the chirality %k = sign(</c| 7^ |^) )- 

The integral over this distribution is normalized to V, 



dZp v x (0 = v. (3.4) 

The index v counts chiralities of the real modes of Dw in Z v : v = % n where n runs over all real 
modes. In the limit of small a the probability of finding configurations with real modes that have 
chiralities of different signs vanishes. In that limit v is simply the number of real modes. The non- 
positivity of the density of real modes is unrelated to this: a change of sign occurs at £ = m. Only 
when m is on the order of or less than 8<3 2 does this have significance in the density since otherwise 
the density is very small anyway. We show an example of the distribution of the chiralities over the 
real modes in fig. 1. 

We now wish to compute the spectrum of the Hermitian Wilson Dirac operator D5 = 75 (Dw + 
m). To that end, introduce the new resolvent 

G v {z,m) = run^lnZj| 1 (m,m,m,0,z,l') = ( Tr f ^q^ jy ( 3 - 5 ) 

and take the discontinuity across the real line. This gives us the spectral density of D5 : 

p?(x) = -Im[G v (jc 5 m)l. (3.6) 
n 

Let us first consider the spectrum corresponding to v = 0. In fig. 2 we show the density 
for fixed m = 5 and various values of a 2 . When a is small, a gap clearly opens up around ±m, 
as it should. The spectrum of D5 then approaches the standard spectrum of the continuum Dirac 



operator with one massive fermion [13], up to a trivial change of variables. In contrast to the 



5 



Microscopic Wilson Dime spectrum for Nf = 1 



Poul H. Damgaard 




-6 -4 -2 2 4 6 



Figure 2: Spectral density of the Hermitian Wilson Dirac operator for v = 0. (Dashed lines Nf = 0.) 



quenched spectrum, the microscopic spectrum of D5 in this Nf =1 theory always has a zero at the 
origin. This is clearly visible in fig. 2. 

There are other differences with the quenched spectrum. Because of the real modes, the spec- 
trum of D5 can change sign in the Nf = 1 theory. A negative density simply corrersponds to a 
theory with a sign problem: the Boltzmann weight in the path integral is not positive definite. The 
existence of a negative density is thus a potential problem for numerical simulations. Fortunately 
the sign problem in this theory is mild: it is only signifact in the small-m limit, and it can be post- 
poned by going to smaller lattice spacings a. We illustrate this phenomenon in fig. 3, where we 
consider a case with v = 1. The analytical understanding we can provide here should be valuable 
for numerical simulations. 

4. Conclusions 

We have presented an explicit computation of the microscopic eigenvalue distributions of the 
Wilson Dirac operator, the real modes of this operator, and the eigenvalues of the Hermitian Wilson 
Dirac operator. We have focused on effects that most clearly distinguish a theory with dynamical 
quarks from the quenched counterpart [p. A sum over the index v can be done straightforwardly. 
This will be presented elsewhere. 
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Figure 3: Same as fig. 2, but now for v = 1. The spectral density is no longer positive. 
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